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ABSTRACT 

Over the past decades open clusters have been the subject of many studies. Such 
studies are crucial considering that the universality of the Initial Mass Function is 
still a subject of current investigations. Praesepe is an interesting open cluster for the 
study of the stellar and substellar mass function (MF), considering its intermediate 
age and its nearby distance. Here we present the results of a wide field, near-infrared 
study of Praesepe using the Data Release 9 (DR9) of the UKIRT Infrared Deep Sky 
Survey (UKIDSS) Galactic Clusters Survey (GCS). We obtained cluster candidates of 
Praesepe based on a 3<r astrometric and 5 band photometric selection. We derived a 
binary frequency for Praesepe of 25.6±3.0% in the 0.2-0.45 M Q mass range, 19.6±3.0% 
for 0.1-0.2 M Q , and 23.2±5.6% for 0.07-0.1 M Q . We also studied the variability of the 
cluster candidates of Praesepe and we conclude that seven objects could be variable. 
We inferred the luminosity function of Praesepe in the Z— and J— bands and derived 
its MF. We observe that our determination of the MF of Praesepe differs from previous 
studies: while previous MFs present an increase from 0.6 to 0.1 M©, our MF shows 
a decrease. We looked at the MF of Praesepe in two different regions of the cluster, 
i.e. within and beyond 1.25°, and we observed that both regions present a MF which 
decrease to lower masses. We compared our results with the Hyades, the Pleiades and 
a Per MF in the mass range of 0.072-0.6 M Q and showed that the Praesepe MF is more 
similar to a Per although they are respectively aged ~85 and ^600 Myr. Even though 
of similar age, the Praesepe remains different than the Hyades, with a decrease in the 
MF of only ~0.2 dex from 0.6 down to 0.1 Mq, compared to ~1 dex for the Hyades. 

Key words: Techniques: photometric - stars: low-mass, brown dwarfs; stars: lu- 
minosity function, mass function - galaxy: open cluster and associations: individual 
(Praesepe) - infrared: stars. 



1 INTRODUCTION 

Over the past decades, open clusters have been the 
subject of many studies (e.g. iBastian et al.l |2010| . and 
references therein). Such studies have brought new 
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sight s into brown dwarf formation (e 
120071 : iBoudreault fc Bailer- Jonesl 12009 
on the discovery of young L and T dwarf and free- 
floating planets (e.g. lBouvier et alJl200sl ; lLodieu et alJl200St 
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Zapatcro Osorio ct al. 2008; Bihain ct al. ll2009l:lLucas et all 
20101 : lOuanz et al.ll2010l : lPena Ramfrez et al.ll2008h . and on 
our understanding of the stellar /subste llar mass function 
(MF) (see review bv lBastian et alj|2010h and their popula- 
tions in the Ga lactic field and in open clusters (see review by 
IChabrierir2 003'). Such studies are crucial considering that the 
universality of the Initial Mass Fu nction (IMF) is still a sub- 
ject of current investiga t ions fe.p.|Kroupa|l2002i;ICovev et aL 
20 111 : iMvers et al1l201ll : iLeigh et alJl2012tlMarks fc Krouua 
20121 ). Most works on the substellar MF have focused on 
young open clusters with ages less than ~150Myr. This is 
partly because brown dwarfs (BDs) are bright when they 
are young, thus aiding the detection of the least massive ob- 
jects. However, the extension of MF studies to older clusters 
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is vital as it allows us to study the intrinsic evolution of BDs 
and how the stellar and substellar population itself evolves. 

Praesepe (M44, NGC2632, RA=8 ,i 40.4 m , 
DEC=+19°4l') is an interesting open cluster to study 
the MF in the stellar and substell ar regimes, c onsid - 
ering its age (r = 590±^ Myr; iFossati et all 120081 ) 
and its distance {(m — M)p = 6.30 ± 0.07 mag, 
d = 181.97±5.7? pc; Ivan Leeuwenl 120091 '). its known 
proper motion (fi a = —35.81 ± 0.29 mas/yr and 
Us — —12.85 ±0.24 mas/vr; lvan Leeuwenl l2009). and the low 
extinction to wards this cluster (E(B - V) = 0.027 ± 0.004; 
iTavlorl I2006T ). while determinations of the metallicity of 
Praesepe go from s o lar-ty pe with [Fe/H]= 0.038±0.039 dex 
l|Friel fc Boesgaardl I1992J - ) to slightly m etal rich with 
[Fe/H]= 0.27±0.10dex (Pace etldl 120081 ). So far, several 
surveys for stellar and substellar objects in the open cluster 



Praesepe have been performed (e.g. Jones fc Cudworthl 



1 19831 ; IHodgkin et alii 19991 ; iBaker et aTj|201Ch . Some surveys 
have the advantage of covering a large a rea and using 
prop er motions, but are rather sha llow (e.g. lHamblv et al.l 
1 19951 ; H raus fc Hillenbrand! l2007h , while other surveys 
used deep photometry, bu t lacked wide areal coverage 



Gonzalez-Garci'a et all 120061 ; iBoudreault et alj l20ld; 



IWang et al.ll201ll ) 

The UKIRT In frared Deep Sky Survey (UKIDSS; 
lLawrence et al.l 120071 ') is a deep large scale inf rared survey 
cond ucted with the wide field camera WFCAM (jCasali et al.l 
l2007h on UKIRT (Mauna Kea, Hawai'i). The survey is 
subdivided into five components: the Large Area Survey, 
the Galactic Clusters Survey (hereafter GCS), the Galactic 
Plane Survey, the Deep Extragalactic Survey, and the Ultra- 
Deep Survey. The GCS aims at covering ~1000 square de- 
grees in 10 star forming regions and open clusters down to K 
= 18.4 mag at two epochs. The main scientific driver of the 
survey is to study the IMF and its dependence with environ- 
ment in the substellar regime using an homogeneous set of 
low-mass stars and brown dwarfs over large areas in several 
regions. The UKIDSS GCS is, therefore, a perfect tool to 
study the open cluster Praesepe, considering the large cover- 
age from the UKIDSS Data Release 9 (DR9) with its relative 
deep photometry spanning from J~10.9mag (i.e. ~fl.7M ) 
down to J~19.3mag (i.e. ~55Mj up ), combined with astro- 
metric information. 

Here we present the results of a wide-field near-infrared 
study of the Praesepe cluster using the DR9 of the UKIDSS 
GCS. The paper is structured as follows. First we present 
the dataset used in our analysis (Section [2}, followed by a 
cross-match with previous surveys (Section [3|. Then we ex- 
tract the new stellar and substellar members in Praesepe 
based on our selection criteria (Section [4}. We discuss the 
level of contamination (Section [5]), the multiplicity of low- 
mass Praesepe members (Section [5J), and the variability of 
our cluster candidates (Section [7]). Finally we derive the lu- 
minosity function (LF) and MF of Praesepe (Section [8} . 



2 THE PHOTOMETRIC AND ASTROMETRIC 
SAMPLE 

The UKIDSS GCS DR9 covers ~ 36 square degrees ob- 
served in five passbands (ZY JHK; iHewett et al.l 120061 ) in 
the Praesepe cluster over a region defined by RA=126-134 



degrees and Dec=16. 5-23.0 degrees. Unfortunately, we are 
missing a region of approximately 0.585 square degrees in 
the central region of Praesepe, delineated by (RA,DEC) 
from (130.37,19.3) to (131.2,19.75) and (130.37,19.75) to 
(130.84,20.2) degrees (Fig.rjJ. 

We selected all good quality point sources in Praesepe 
detected in JHK with no constraint on detection in Z, Y , 
and K2 (second epoch K band). We imposed a selection 
on point sources only in JHK and pushed the complete- 
ness towards the faint end by insisting on the ClassStat 
parameters which represent the roundness of an image 
(i.e. — 3<ClassStat<-|-3). The Structured Query Language 
(SQL) query used to select sources along the line of sight of 
the Praesepe i s ident ical to the query used for the Pleiades 
IjLodieu et al.l l2012al). The SQL query includes the cross- 
match with 2MASS (cutri et al.ll2003l ; lskrutskie et al.f2006l ) 
to compute proper motions for all sources brighter than the 
2MASS 5cr completeness limit at J = 15.8 mag as well as 
the selection of proper motion information from that newly 
available in UKIDSS DR9. We used the GCS proper mo- 
tion measurements in this work as they are more accurate 
due to the homogeneous coverage, completeness, and spatial 
resolution of the GCS. We limited our selection to sources 
fainter than Z = 11.7, Y = 11.8, J = 10.9, H = 11.4, Kl 
= 10.0 and K2 = 10.4 mag to avoid saturation. The com- 
pleteness limits, taken as the magnitude where the straight 
line fitting the shape of the number of sources as a function 
of magnitudes falls off, are Z = 20.4, Y = 20.0, J = 19.3, 
H = 18.7, Kl = 17.9, and K2 = 18.4 mag. (The method 
for de riving the detection limit is presented in lLodieu et al.l 
12009.1 ) 

The query returned 218,141 sources with J = 10.9-21.2 
mag over 36 square degrees towards Praesepe. The full GCS 
DR9 coverage is displayed in Fig.[T] The resulting (Z — J,Z) 
colour-magnitude diagram (CMD) is shown in Fig. [2] 

Proper motion measurements are available in the WF- 
CAM Science Archive for UKIDSS data releases from the 
Data Release 9 (DR9) for all the wide/shallow survey with 
multiple epoch coverage in each field, including the GCS. 
The details of the procedure to derive the relative proper 
motions is described in detail in ICollins fc Hamblvl (|2012h 
and a summary is given in lLodieu et al.l ( 2012a '). The proper 
motion vector point diagram for all sources detected towards 
Praesepe is presented in Fig. [3] 



3 CROSS-MATCH WITH PREVIOUS 
SURVEYS 

We compiled a list of member cand idates of Preasepe 
publis he d from previo us stu dies in iJones fc CudwortI 



( 19831), lHamblv et all 



iGonzalez-Garcfa et al 



IHodgkin et alj 



Kraus & Hillen brand! 



j|2007l). IBaker et al.1 (|2010l ). IBoudreault et alj j|2010l ) and 
IWang et al. ( 2011 ). We matched these previously known 
members with the 218,141 sources from UKIDSS GCS DR9 
towards Praesepe. Only those members with positions that 
are within 3 arcsec of the location of each object in these 
previous catalogues are included in our sample. 

In Table Q] we present the numbers of cluster candidates 
published by each study (All) and the corresponding num- 
bers of cluster candidates retrieved by the GCS (DR9). We 
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Figure 1. All detected objects towards Praesepe as released by 
the UKIDSS GCS DR9 (small black dots). For clarity, only one 
of ten objects are plotted. We also show the objects which are 
recovered in our astrometric and photometric selection (larger 
black dots, see Section |4}. The missing region at the centre is due 
to poor quality images that were discarded during the quality 
control process of UKIDSS. 




Z - J [mag] 

Figure 2. CMD of Z vs. Z — J of all detected objects towards 
Praesepe from the UKIDSS GCS DR9. The cluster sequence of 
Praesepe is clearly visible between [Z — J,Z)~(0.6,12.5) and (Z — 
J,Z)~(1.3,17). 
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Figure 3. Proper motion diagram of all detected objects towards 
Praesepe from the UKIDSS GCS DR9. The position of Praesepe 
is clearly visible at fi a = —34.17 ± 2.74 mas/yr and n$ = — 7.36± 
4.17 mas/yr. 



also present in this table the number of objects recovered in 
our astrometric and photometri c selections (see Sect ion H}. 

For some wo rks (e.g. Ijones fc Cudworthl 1 19831 : 
iHodekin et~aT1ll999l l. the ratios of Praesepe member candi- 
dates published in the literature and recovered by the GCS 
DR9 are low (i.e. 51.5% and 58.8% respectively). These 
older works had shallow photometry and most of their 
candidates are saturated in the GCS DR9. Other works 
have large coverage beyond the area of the G CS DR9, with 
more than 30 and 300 square d e grees for iHamblv et al.l 
(1995) and iKraus fc Hillenbrand] i|2007l ). respectively. As 
for the cluster candidates from previous deep surveys 
(i.e. iGonzalez-Garci'a et alJ l200d : iBoudreault et all l2010l : 
IWang et al]|201ll ). the objects not recovered are either too 
faint for UKIDSS GCS, or they are in the central area not 
covered (see Fig. [lj. Finally, all but one (not detected in 
the second epoch of the AT-band observations) object from 
iBaker et ail fcOlCh are recovered, since this work was based 
on an earlier release of the UKIDSS GCS. 

In Fig.|4]we present in four CMDs all Praesepe member 
candidates published in the literature and recovered by the 
GCS DR9. We use these previous cluster candidates with 
proper motions consistent with Praesepe membership (and 
those classified as proper motion non members) to define se- 
lection lines which will be used for our photometric selection 
(see Section Q. 



4 NEW STELLAR AND SUBSTELLAR 
MEMBERS IN THE PRAESEPE 

The selection used in this work and presented in Sec- 
tion [?7T] and 22] JJLjdentical to o ur studies of the Pl eiades 
jLodieu et al]|2012al ) and a Per (|Lodieu et al]|2012bl ) clus- 
ters. The main scientific goals of these studies is to perform 
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Figure 4. CMDs showin g prev i ously reported membe rs of Praesepe r ecovered in UKIDSS DR9 GCS. Previous known members are taken 
fromljones fc Cudwortbl [U983ri.lHamblv et all dl995ll.lHodgkin et all l|l999l l. lGo~nzalez-Garci'a et all l|200cJ l. lKraus fc Hillenbrand! d2007t ). 
Bake r et all i2010h . iBoudreault et all d201(J ) and lWang et all ll201ll) . Filled dots are objects which have a proper motion consistent with 
Praesepe membership (3cr astrometric selection), while objects presented with crosses don't. The NextGen isochrones are the solid blue 
lines with masses on the right, while the DUSTY isochrones are dotted green lines with masses on the left. The dashed red lines represent 
our photometric selection criteria in each diagram as detailed in Section [4] 



an homogeneous analysis of the LF and MF of the Pleiades, 
a Per, and Praesepe. 

4.1 Astrometric selection 

First we performed a selection based on the proper motion o f 
Praesepe. The proper motion values of Ivan Leeuwenl (120091 ) 
are the absolute proper motion measurements obtained with 
Hipparcos. Here, the proper motion values from UKIDSS are 



relative to a reference frame defined by the faintest, most nu- 
merous objects detected in each frame. The position of the 
cluster in the relative proper motion diagram of Fig. [3] is 
^ Q = -34.17 ± 2.74 mas/yr and fj, s = -7.36 ± 4.17 mas/yfl 
We select all objects within 3 a^ Q of \x a and within 3 a^ s of 

1 We simply used the median of the distribution in fi a (between 
—40 and —30 mas/yr) and fis (between —15 and —2 mas/yr) for 
all objects within J<17mag. 
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Table 1. Updated membership of the Praesepe member 

candidates published in the literature and recovered by 

the GCS DR9. Papers dedicated to Praesepe over the 

past two decades are ordered by year. References are: 



Jones & Cudworth (1983), Hamblv ct al. (1995), Hodgkin ct al. 


(1999 


), 


Gonzalez-Garcfa et al.l d2006f). iKraus & Hillenbrand 


(2007 


Baker et al.l ||2010|). iBoudreault et alj d2010h and 



IWang et al] l|201ll ). Column 2 gi 



gives the numbers of cluster 
member candidates published in the reference given in column 
1. Column 3 gives the numbers of sources recovered in GCS 
DR9. Column 4 gives the percentages of sources recovered in the 
GCS DR9. Column 5 and 6 give the number of members and 
non-members based on our photometric and astrometric analysis 
respectively. 



Survey 


All 


DR9 


% 


Memb 


NM 


Jones & Cudworth (1983) 


206 


106 


66.5 


49 


57 


Hamblv et al. (1995) 


515 


383 


74.4 


323 


60 


Hodgkin et al. (1999) 


17 


10 


58.8 


6 


4 


Gonzalez-Garcfa et al. (2006) 


20 


11 


70.0 


1 


13 


Kraus & Hillenbrand (2007) 


1130 


711 


62.9 


654 


57 


Baker et al. (2010) 


147 


146 


99.3 


125 


21 


Boudreault et al. (2010) 


150 


100 


66.7 


37 


61 


Wang et al. (2011) 


59 


29 


49.2 


7 


22 



Us , where a l _ la and 3 <r MJ are the error on the measurement of 
Ha and [is of each source in the UKIDSS GCS DR9 respec- 
tively. All the objects selected astrometrically within 3 a are 
presented in a CMD of Z vs. Z - J in Fig. From 218,141 
sources detected in the UKIDSS GCS DR9 towards Prae- 
sepe, 33,854 objects are retrieved from our 3u astrometric 
selection (15.5%). 



4.2 Photometric selection 

We used the selection lines plotted in the CMDs of Fig.|4]to 
perform the photometric selection of our cluster candidates. 
These lines were chosen based on the position of previously 
reported Praesepe members recovered in UKIDSS GCS DR9 
and by our astrometric selection. These lines are ~2.5 sigma 
bluer than the location of the cluster sequence, where the 
value of sigma is the distribution in colour of the cluster se- 
quence in four different bins of magnitude. The photometric 
criteria were defined as follows: 

• Candidates should be detected in the ZYJHK- 
bands 

• Candidates should lie above the lines defined by (Z — 
J,Z) = (0.4,11) and (1.1,17.5), and by (Z - J ,Z) = (1.1,17.5) 
and (2.5,21.5) 

• Candidates should lie above the lines defined 
by (Z - K,Z) = (0.8,11) and (1.85,18.5), and by (Z - 
K ,Z) = (1.85,18.5) and (3.7,21.75) 

• Candidates should lie above the lines defined 
by (Y - J,Y) = (0.3,11.5) and (0.55,17.5), and by (Y - 
J,Y) = (0.55,17.5) and (1.45,21) 

• Candidates should lie above the lines defined 
by (J - K,J) = (0.65,11) and (0.75,18.25), and by (J - 
K,J) = (0.75,18.25) and (1.65,20) 

This is presented in Fig. [6] where we show the same four 
CMDs in Fig.Uand the objects recovered by our astrometric 
selection. We also present the previously reported members 
of Praesepe recovered in UKIDSS DR9 GCS. 
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Figure 5. CMDs of Z vs. Z — J of all detected objects towards 
Praesepe from the UKIDSS GCS DR9 with a proper motion 
within 3<t of the cluster's mean proper motion. The masses cor- 
responding to the Z magnitudes are given on the right (using a 
NextGen isochrone) and on the left (using a DUSTY isochrone). 
The black and white dashed lines are our photometric selection 
criteria (see Fig. |4") . 



Table 3. Praesepe member candidates based on our astrometric 
and photometric selection. Here we present the physical parame- 
ters of the candidates, i.e. mass and effective temperature (T c ff). 
The nomenclature used to identify the objects in this table is 
UGCS JHHMMSS.SSiDDMMDD.D. Therefore, the first object 
in the table is UGCS J084644.83+221246.3. Table[3]is published 
in its entirety in the electronic edition of the MNRAS. A fraction 
is shown here for guidance regarding its form and content. 



RA 

/h m s\ 



DEC 

(° ' ") 



Mass 
[M©] 



T c ff 
[K] 



8:46:44.83 +22:12:46.3 0.132±0.005 3062±125 

8:43:15.16 +22:01:47.2 0.198+0.013 3286+217 

8:34:04.31 +16:58:24.7 0.269+0.006 3397+ 74 

8:37:31.99 +16:57:35.7 0.336+0.013 3465+138 



From 33,854 objects from our astrometric selection, 
1,116 are retrieved from our photometric selection (3.3%). 
This final selection gives us the objects which will be used 
for our following analysis. Our cluster candidates are listed 
in Tableland in Table [3] 



4.3 Search for lower mass members 

Here we search for fainter and cool substellar members of 
the Praesepe by trawling our selection first for stars having 
no Z-band detection, and secondly for stars having no Z 
and Y-band detections. 
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Figure 6. Same as Fig. [4] showing previously reported members of Praesepe (small dots) recovered in UKIDSS DR9 GCS, and for which 
the proper motion is consistent with Praesepe membership. The larger dots are objects for the UKIDSS DR9 GCS, but not reported 
before in the literature as cluster members. From one panel to the other, we only keep the objects redder than the red dash lines, starting 
from Z vs. Z — J (top-left) followed by Z vs. Z — K , to Y vs. Y — J and finally with J vs. J — K. 



4-3.1 Search for lower mass members detected in YJHK 

To extend the cluster sequence to fainter objects, we 
searched for Praesepe candidate members undetected 
in Z. We imposed similar photometric and astrometric 
criteria as in Section r4.1l and l4~2l but without detection in Z: 

• Candidates should not be detected in the Z-band 

• Candidates should have a proper motion within 3a of 
the cluster's mean proper motion 

• Candidates should have Y^18 and J^18.8mag 



• Candidates should lie above the lines defined by 
(Y - J,Y) = (0.3,11.5) and (0.55,17.5), and by (Y - 
J,Y) = (0.55,17.5) and (1.45,21) 

• Candidates should lie above the lines defined by 
(J - K,J) = (0.65,11) and (0.75,18.25), and by (J - 
K,J) — (0.75,18.25) and (1.65,20) 

Out of 17 objects selected based on astrometry and 
without detections in Z, 16 are selected based on the colour 
cut above (94.2%). All 16 objects have an existing image in 
the Z-band in the UKIDSS GCS DR9 database. A list of 
these 16 objects is presented in Table [4] 
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Table 2. Praesepe member candidates based on our astrometric and photometric selection. Here we present the photometry and the 
astrometry of the candidates. The nomenclature used to identify the objects in this table is UGCS JHHMMSS.SS+DDMMDD.D. 
Therefore, the first object in the table is UGCS J084644. 83+221246. 3. Tablc[2]is published in its entirety in the electronic edition of the 
MNRAS. A fraction is shown here for guidance regarding its form and content. 



RA 

/h m s\ 



DEC 

(° "') 



Z 
[mag] 



Y 
[mag] 



J 

[mas 



H 

[mag] 



Kl 

[mag] 



K2 
[mag] 



fla COS 8 

[mas/yr] 



[mas / yr] 



1:46:44.83 +22:12:46.3 16.914+0.011 16.449+0.009 15.865+0.009 15.310+0.009 15.024+0.012 15.029+0.006 -26.85+ 3.54 -16.99+ 3.54 
1:43:15.16 +22:01:47.2 16.005+0.006 15.546+0.005 15.024+0.005 14.427+0.005 14.143+0.005 14.162+0.003 -29.28+ 3.55 -3.86+ 3.55 



1:34:04.31 
1:37:31.99 



-16:58:24.7 
-16:57:35.7 



15.394+0.005 
14.978+0.004 



14.889+0.004 
14.522+0.003 



14.321+0.004 
13.905+0.003 



13.757+0.003 
13.356+0.002 



13.455+0.003 
13.052+0.002 



13.454+0.002 -31.56+ 3.45 2.76+ 3.45 
13.051+0.002 -34.85+ 3.45 -3.12+ 3.45 



Table 4. Praesepe member candidates based on our astrometric and photometric selection using Y JHK photometry (top lines) and using 
JHK photometry (last two lines). The nomenclature used to identify the objects in this table is UGCS JHHMMSS.SS+DDMMDD.D. 
Therefore, the first object in the table is UGCS J084322. 72+222646.6. 



RA 

Ah m e 


) 


DEC 

(° ' ") 




Z 
[mag] 




Y 
[mag] 






J 

[mag] 






H 

[mag] 






Kl 
[mag] 






K2 
[mag] 




/Ma COS 8 

[mas/yr 




US 
[mas / yr] 


8 


43 


22 


72 


+22:26 


46 


6 




19 


734+0 


123 


18 


154+0 


049 


17 


377+0 


062 


16 


630+0 


049 


16 


641+0 


029 


-34.55+ 5 


91 


-16.26+ 5.91 


8 


45 


07 


82 


+21:17 


29 


2 




19 


665+0 


146 


18 


358+0 


082 


17 


770+0 


063 


17 


041+0 


058 


17 


041+0 


041 


-11.85+ 7 


63 


-1.27+ 7.63 


8 


40 


20 


81 


+20:05 


21 


4 




19 


310+0 


076 


18 


241+0 


049 


17 


637+0 


072 


16 


966+0 


072 


16 


898+0 


033 


-43.70+ 5 


78 


-9.07+ 5.78 


8 


37 


15 


19 


+20:09 


51 


2 




19 


314+0 


076 


18 


097+0 


047 


17 


293+0 


039 


16 


553+0 


049 


16 


560+0 


025 


-34.20+ 5 


97 


1.55+ 5.97 


8 


36 


36 


89 


+20:07 


00 


7 




19 


269+0 


072 


18 


225+0 


053 


17 


433+0 


044 


16 


717+0 


054 


16 


727+0 


029 


-29.03+ 6 


38 


-8.23+ 6.38 


8 


38 


03 


30 


+19:47 


15 







19 


881+0 


147 


18 


333+0 


068 


17 


403+0 


043 


16 


720+0 


049 


16 


725+0 


031 


-31.72+ 8 


39 


7.77+ 8.39 


8 


45 


■17 


32 


+20:05 


25 


8 




19 


686+0 


117 


18 


504+0 


076 


17 


848+0 


069 


17 


534+0 


110 


17 


570+0 


068 


-13.04+ 8 


35 


-27.58+ 8.35 


8 


34 


22 


13 


+19:54 


38 


7 




20 


011+0 


161 


18 


623+0 


087 


18 


104+0 


094 


17 


560+0 


124 


17 


536+0 


068 


-28.77+10 


31 


-34.87+10.31 


8 


39 


22 


81 


+19:40 


37 


1 




19 


743+0 


113 


18 


472+0 


066 


17 


751+0 


071 


17 


050+0 


068 


17 


123+0 


039 


-42.70+ 6 


02 


-2.48+ 6.02 


8 


39 


01 


37 


+20:35 


04 


1 




19 


735+0 


111 


18 


546+0 


053 


17 


882+0 


071 


17 


286+0 


070 


17 


280+0 


051 


-38.44+ 6 


47 


-4.31+ 6.47 


8 


38 


05 


45 


+19:23 


27 


7 




19 


484+0 


095 


18 


246+0 


057 


17 


510+0 


045 


16 


993+0 


057 


16 


882+0 


032 


-39.10+ 8 


34 


-14.61+ 8.34 


8 


17 


35 


00 


+17:49 


53 







19 


925+0 


149 


18 


602+0 


099 


17 


854+0 


065 


17 


409+0 


085 


17 


523+0 


060 


-21.59+11 


13 


-14.85+11.13 


8 


44 


25 


06 


+17:50 


52 


9 




19 


743+0 


124 


18 


418+0 


085 


17 


799+0 


063 


16 


883+0 


052 


16 


930+0 


035 


-28.17+ 6 


92 


8.40+ 6.92 


8 


54 


12 


81 


+18:26 


11 


7 




19 


825+0 


146 


18 


567+0 


075 


17 


943+0 


072 


17 


377+0 


099 


17 


300+0 


050 


-30.57+ 8 


63 


-2.17+ 8.63 


8 


35 


32 


57 


+17:06 


30 


9 




19 


562+0 


105 


18 


436+0 


064 


17 


852+0 


073 


17 


193+0 


071 


17 


160+0 


047 


-14.73+ 6 


72 


-18.78+ 6.72 


8 


33 


42 


21 


+18:47 


10 


3 




19 


784+0 


110 


18 


427+0 


066 


17 


615+0 


061 


17 


015+0 


064 


17 


029+0 


039 


-33.11+ 6 


29 


-9.67+ 6.29 



8:47:13.71 +18:16:38.3 18.693+0.119 17.675+0.060 16.954+0.064 16.967+0.041 -40.67+ 9.93 6.83+ 9.93 

8:43:52.01" +17:39:20.7° 18.552+0.091 17.894+0.066 17.178+0.063 17.203+0.044 -19.91+ 9.48 -17.76+ 9.48 



a A visual inspection of the ZY J H K1K2 images suggests a possible companion. 



4-3.2 Search for lower mass members detected in JHK 



4.4 Formal Membership probabilities 



We repeated the procedure described above looking for Z 
and Y non detections. We applied the following criteria: 



• Candidates should not be detected in Z and Y-bands. 

• Candidates should have a proper motion within 3a of 
the cluster's mean proper motion 

• Candidates should have J^18 and J^18.8mag. 

• Candidates should lie above the lines defined by 
(J - K, J) = (0.65,11) and (0.75,18.25), and by (J - 
K,J) = (0.75,18.25) and (1.65,20). 

Out of 71 objects selected based on astrometry and 
without detections in Z and Y, 17 are selected based on 
the colour cut above (23.9%), but only 2 objects have an 
existing image in the Z and Y-band in the UKIDSS CCS 
DR9 database and remain undetected. These 2 objects are 
listed in Table H 



As a check on our photometric and astrometric selection, we 
calcul ate formal membership pro babilities based on the anal- 
ysisof Deacon fc HambTvl (|2004l ) . This uses a method similar 
to lSandersI l|l97ll ) who first proposed calcul ating member- 
ship p robabilities from proper motions and to lHamblv et al.l 
(|l995l ) who built on this method. In this method the proper 
motion vector point diagram is rotated so that the cluster 
lies on the y' axis, 



x — cos(— 0.557r)a; — sin(— 0.557r)y 
y — cos(— 0.557r)j/ + sin(— 0.557t)j; 



(1) 
(2) 
(3) 



and the distribution of proper motions in this vector point 
diagram is defined as 

* = + (1 - /)#c (4) 
where is the field star distribution defined as 



Co 



V / (27r)E ; 



■<**(.- f 2g ) 



(5) 
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Put simply the field population is modelled as an expo- 
nential in the y direction with scale length r, a gaussian in 
the x' direction with mean fi x f and standard deviation E x 
and has normalisation cq. $ c is the cluster distribution 



■exp- 



2^2 



(6) 



The cluster is fitted with a bivariant gaussian with a single 
standard deviation (a) centred on fi xc and fj, yc . The mem- 
bership probability is defined as, 

m 

In iDeacon fe Hamblvl l|2004l ) and lLodieu et aD l|2007bl ) we 
allowed all our parameters to be freely fitted by a maxi- 
mum lik elihood method. How ever here we take the approach 
taken in lLodieu et all (|2012al ) where the standard deviation 
of the cluster proper motions is fixed to be the proper motion 
measurement errors calculated in the appropriate magnitude 
bin. 

The input list for the fitting procedure had the same 
source quality parameter cuts as listed in Section [5] with a 
single colour cut being applied, 

• Z < 7 x (Z - J) + 9.5 & Z < 18.6 

• Z <3x(Z - J) + 14&Z > 18.6 

The fitted parameters for each brightness range are 
shown in Table [5] Note that the cluster position appears to 
migrate in the last fitted bin. This is likely not due to any 
real variation in the cluster proper motion but to the fitting 
procedure becoming unreliable for the faintest objects. 

Table[5]contains a truncated list of high priority (p>0.5) 
candidates, the full list is available in electronic form. Us- 
ing the probabilistic analysis, we obtained a total of 1,015 
candidates with a membership probability higher than 50%, 
comparable to the 1,116 cluster candidates based on a pho- 
tometric and astrometric selection. These objects are shown 
in Fig. [5] (right panel) in the CMD of K vs. Z - K with all 
GCS point sources towards Praesepe. 



5 ESTIMATION OF CONTAMINATION 

The number density of field objects in our final list of can- 
didat es as a functi o n of m ass was obtained in a similar way 
as m IJeffries et all (|2004r i. We obtained the radial profile of 
our cluster candidates in three mass ranges: above 0.4 Mq, 
between 0.15 and 0.4 M Q , and below 0.15 M (Fig. [7). 

However, we must consider the first and last two dat- 
apoints with caution, as we are not complete in coverage 
within 1° and beyond 3°. Therefore, we are not complete 
at the tidal radius of Prae sepe at 3.5°±0.1° (11.5±0.3pc; 
iKraus fe Hillenbrand! [200^}_ and not within its core radius 
at 1.1°±0.1° (3.5±0.3pcT lAdams et al]|2002h . 

Therefore, we only use the bins between 2.5 and 3deg. 
at each mass range to obtain an upper limit of contami- 
nation. This gives 1.27 objects per square degree for can- 
didates with masses above 0.4 Mq, 3.70 between 0.15 and 
0.4 Mq, and 2.08 objects per square degree for candidates 
with masses below 0.15 Mq. Within 3deg. from the clus- 
ter centre, this gives a contamination of 11.9%, 9.8% and 
23.8,% for the same mass range respectively, or 18.7% for 
the whole Praesepe sample. 

These numbers appear quite large. We stress again that 
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Figure 7. Radial density plots of our candidate members of Prae- 
sepe in three mass ranges: above 0.4 Mq (top panel), between 0.15 
and 0.4 Mq (middle panel), and below 0.15 Mq (low panel). The 
error bars on each datapoint arc Poissonian arising from the num- 
ber of objects in each bin. The arrows represent lower limits, as 
we are not complete in these bins. 



these are upper limits, since the coverage is not complete 
at the tidal radius of Praesepe and beyond. However, we 
can claim completeness in our cluster candidate list, and 
MF determination, better than 90 %. This is justified by the 
fact that our astrometric selection includes all objects within 
3cr of the cluster's mean proper motion (completeness of 
>99%) and that the lines used in our photometric selection 
go at least 2a bluer from the cluster main sequence in all 
the CMD used for the photometric selection (completeness 
of ~95.4%). 



© 2012 RAS, MNRAS 000,[THl7] 



IMF from UKIDSS: III Praesepe 9 



Table 5. A table containing the fitted parameters for each magnitude range. For definitions of each term see lDeacon fc Hambl 

Note in this case we did not freely fit the cluster proper motion standard deviation a. Instead this was calculated from the mean proper 

motion error in each magnitude bin. 



Zrange N sta rs f cr fi xc fj, yc t T, x fi x f 



12-13 


272 


0.821 


3.42 


-2.03 


35.63 


13.97 


17.01 


-4.60 


13-14 


537 


0.743 


3.40 


-2.55 


35.24 


16.67 


16.95 


-7.73 


14-15 


866 


0.778 


3.41 


-2.31 


35.04 


14.89 


16.64 


-6.41 


15-16 


1146 


0.823 


3.40 


-2.31 


35.00 


13.78 


16.66 


-6.69 


16-17 


1251 


0.844 


3.48 


-2.79 


35.02 


11.46 


15.42 


-6.81 


17-18 


906 


0.936 


3.67 


-2.34 


35.08 


11.24 


14.11 


-6.50 


18-19 


587 


0.936 


4.37 


-3.75 


34.93 


11.44 


14.33 


-5.96 


19-20 


316 


0.850 


6.04 


-2.51 


34.16 


11.88 


16.05 


-6.19 


20-21 


241 


0.863 


9.02 


-3.17 


32.83 


16.02 


18.09 


-4.58 



Table 6. Praesepe member candidates identified to have membership probabilities calculated from proper motions greater than 0.5. The 
nomenclature, as used above to identify the objects in this table, is UGCS JHHMMSS.SSiDDMMDD.D. Table [6] is published in its 
entirety in the electronic edition of the MNRAS. A fraction is shown here for guidance regarding its form and content. 



Name 


Membership 
probability 


Z 


Y 


J 


H 


K 


UGCS J084344.72+211234.2 
UGCS J084037.87+202017.9 


0.85 
0.92 


12.996±0.001 
12.865±0.001 


12.710±0.001 
12.614±0.001 


12.203±0.001 
12.169±0.001 


11.615±0.001 
11.607±0.001 


11.360±0.001 
11.379±0.001 


UGCS J083814.20+194723.4 
UGCS J083829. 61+195145.0 


0.83 
0.82 


12.464±0.002 
12.195±0.001 


12.189±0.001 
12.010±0.001 


11.744±0.001 
11.695±0.001 


11.500±0.001 
11.617±0.001 


10.942±0.001 
10.918±0.001 



In Fig.[8](left panel) we present the K vs. Z-K CMD of 
all GCS point sources towards Praesepe (small dots) and the 
cluster candidates based on our astrometric and photomet- 
ric selections (larger dots). We observe three structures in 
this diagram. The structure at Z — 7f~0.6mag contains pre- 
dominantly Galactic disk turn-off stars and the structure at 
Z — A"~1.7 mag represents Galactic disk late- type and giant 
stars, while the structure at Z — K<^2.5 mag is mostly com- 
posed of galaxies. This indicates that most of the contami- 
nants of our cluster candidates with masses above 0.15 Mq 
would be Galactic disk late-type and giant stars, while most 
of the contaminants of candidates less massive than 0.15 Mq 
would include Galactic disk late-type and giant stars, but 
also unresolved galaxies. 



6 DISCUSSION ON THE BINARY 
FREQUENCY 

Here we investigate the multiplicity of low-mass and very 
low-mass Praesepe members using the photometry and 
colours from the GCS. 

We selected binary candidates in the following way. 
We used the Z vs. Z — K CMD (Fig. 0) to define three 
mass ranges: between 0.45 and 0.2 Mq, between 0.2 and 
0.1 Mq, and between 0.1 and O.O7M0. These mass cuts 
were chosen as they can be compared directly with the 
binary fracti o n obt ained from the hydrodynamical simula- 
tion bv lBatel (|2012T l (with the mass intervals of 0.5-0.2 M Q , 
0.2-0.1 M Q and 0.1-0.07 M ). At each transition of these 
mass ranges, we defined an area starting from 0.376 mag 
brighter to 1.13 mag fainter than the cluster sequences (sin- 
gle objects), and an area starting from 0.376 mag brighter 



to 1.5 mag brighter than the cluster main sequence (binary 
candidates). The binary fraction is the ratio between the 
number of objects in the brighter area, with the sum of both 
areas, for each mass range. Because of our 3<r selection in 
astrometry and ~2.5<r bluer in photometry, we expect to be 
complete up to 97% in our selection between the saturation 
and detection limit, bracketing the mass intervals used for 
the binary fraction estimation. 

Our binary candidates are presented in a Z vs. Z — K 
CMD shown in Fig. [5] The binary fraction we obtain is 
25.8±3.7% in the mass range 0.07 to 0.1 M Q , 19.6±3.8% in 
0.1 to 0.2 Mq, and 23.3±7.3% in 0.2 to 0.5 M . We present 
these binary fractions in Table and we compare them di- 
rectly wit h tho s e obt ained from the hydrodynamical simu- 
lation by Bate (2012] 2 ], and a lso the photometric determi- 
nation bv lPinfield et al l |2003h . 

We observe a similar binary fraction to lBatd (|2012t ) be- 
tween 0.07-0.1 Mq (within la, where a is the error on our 
measurements), a disagreement of only ~1.5<r for the binary 
fraction in the mass range 0.2-0.45 Mq, but a significantly 
lower binary fraction by more than 3<r for the mass range 
0.1-0.2 M . However , we note that the binary fraction in 
Fig. 7 of iBatd |2012l ) is also higher compared to previous 
observat ions of the binary fraction for field dw arfs below 
0.2 Mq (|Basri fc Reiners!l2006l ; Iciose et al.ll2003T l. 

Our binary f ractio n is lower by ~1.8<t than the value 
of lPinfield ei~al] (|2003t ) for masses in the range of 0.35 to 
0.6 M Q (0.2 to 0.45 M for our work). On the other hand, 

2 Table 2 from (Bate 2012) presents the number of single objects, 
and also the number of binary, triple and quadruple systems per 
each mass range obtained from their hydrodynamical simulation. 
In our analysis, we only consider the binary system. 
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Figure 8. Left panel: CMD of K vs. Z — K of all detected objects towards the Praesepe from the UKIDSS GCS DR9. For clarity, 
only one of five objects are plotted. We also show the objects which are recovered in our astrometric and photometric selection (larger 
green dots, see Section [!}• The structure at Z — i\~0.6mag contains predominantly Galactic disk turn-off stars and the structure at 
Z — K~1.7mag is Galactic disk late-type and giant stars, while the structure at Z — A"~2.5mag is mostly composed of galaxies. The 
NextGen isochrone masses corresponding to the K magnitudes are on the right, while the DUSTY isochrone masses are on the left. Right 
panel: Same as the left panel, but using the objects with a membership probability higher than 50% from the probabilistic analysis. 



their binary fraction for lower masses is significantly higher 
by a factor of two than our values (more than 4<r for the mass 
range 0.07 to 0.1 Mq, and more than 8<r for the mass range 
0. 1 to . 2 Mm) and are also significantly hig her than those 
of Bate! |2012h (considering the error bars of Pinfiel d et al.l 
2003). We attribute this overes timation of binary fre quency 
to the selection method used bv lPinfield et~aH (|2003r i. These 
authors used a method similar to the one used here, i.e. by 
defining an area above the cluster sequences corresponding 
to the binary population. However, they used a cluster se- 
quence constructed both from isochrones and empirical data. 
The resulting assumed cluster sequence appears bluer than 
the observed cluster sequ ence by about I — K ^0.2 mag at 
if~14mag (see Fig. 3 of IPinfield et all 120031 '). This would 
give a fainter equal-mass binary sequence, and therefore, 
overestimate the number of binaries and the binary fraction. 

Using the probabilistic analysis presented in Section l4~4l 
we obtain similar binary fractions using objects with a mem- 
bership probability higher than 50%, i.e. 27.9% in the mass 
range 0.07 to 0.1 M , 18.0% in 0.1 to 0.2 M Q , 26.6% in 0.2 
to 0.5 M©, and 23.7% in the whole mass range from 0.07 to 
O.5M . 



7 VARIABILITY AT 600 MYR 

Here we discuss the variability of the Praesepe member can- 
didates using the two epochs provided by the GCS in the 
K-band (Kl being the first epoch, K2 the second). 

In Fig. [TO] we present the difference in the K magnitudes 
(i.e. K1-K2) as a function of the K magnitude (average of 
Kl and K2) for all the photometric and astrometric can- 



Table 7. Binary fracti o n from our work on Praesepe compared to 
the simulation of lBatel J2012T) . Errors on the fractions are Poisso- 
nian arising from the number of single objects and the number of 
binary candid a tes. W e also give the binary fraction observed by 
IPinfield et al. I (I2003T) . Numbers obtained using the probabilistic 
analysis are in parenthesis. 



Mass range (Mq) 


Single 


Binary 


% 


Bate 


0.07 - 0.1 
0.1 - 0.2 
0.2 - 0.45 


73(44) 
225(201) 
279(256) 


22(17) 
55(44) 
96(91) 


23. 2±5. 6(27.9) 
19.6±3.0(18.0) 
25.6±3.0(26.2) 


27.3±11.6 
29.2±5.3 
32.1±4.7 a 


0.07 - 0.45 


588(538) 


176(167) 


23.1±2.0(23.3) 


29.2±2.0 


Pinfield et al. 
0.09 - 0.11 
0.2 - 0.35 
0.35 - 0.6 






47 +i3 
44±6 
31±S 





a I Bate! J2012f) used the interval 0.2 to 0.5 M , but the values 
presented here are for the interval 0.2 to 0.45 M (Bate 2012, 
private communication). 

didates. The brightening observed at _R"<12mag is due to 
the difference in depth between the first and second epoch, 
which causes a decrease of 0.5 mag in the completeness limit, 
but also in the saturation limno Therefore, objects with 
if<12mag are excluded from our variability analysis. 

We selected as variable all objects with \K 1-K2\ larger 



3 The exposure times have been doubled for the second epoch 
with relaxed constraints on the seeing requirement and weather 
conditions. 
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Figure 9. CMD of Z vs. Z — K of all our cluster candidates of Praesepe from the UKIDSS GCS DR9. The NextGen isochrone masses 
are on the right, while the DUSTY isochrone masses are on the left. Left panel: The red dotted boxes show the selection areas for single 
objects (bottom) and for binaries (top). Right panel: Dots are single objects and circles with dots are binary candidates with masses 
between 0.5 and 0.07 Mq. 



than 3a, where a is the standard deviation of KX-K2 around 
the K magnitude of the object. Nine objects are identified as 
variable candidates. We looked at the UKIDSS GCS DR9 
ZYJHK1K2 images and give our comments in Table [8] 
One source has its photometry affected by the presence of a 
nearby object, while one other might show spurious variabil- 
ity due to poorer image quality in one of the isT-band images. 
There are seven remaining objects which are likely to be vari- 
able objects in Praesepe, which is 0.63% of our member can- 
didates (black circles with dots in Fig. I10[) . One of the vari- 
able candidates, UGCS J084353.41+210126.3, is in the sub- 
stellar regime with an estimated mass of 0.064T0.013 Mq. 

Using the probabilistic analysis, we also obtain 6 objects 
with \K\-K2\ larger than 3er (0.59% of the cluster candidate 
having a membership probability higher than 50%, including 
3 from the variable candidates sample from the astrometric 
and photometric selection). 



8 LUMINOSITY AND MASS FUNCTION 
8.1 Isochrone of Praesepe 

We derived a NextGen isochrone u sing a grainless at- 
mosphere and evolut ionary tracks (Ba raffe et al.l 1 19981 ; 
lHauschildt et al.l Il999f l and a DUSTY iso chrone using a 
dusty atmosphere and e volutionary tracks |Chabrier et al.l 
|2000| ; lAllard et ai"1l200li ). We used the known astrophysical 
parameters of Praesepe (age, distance) presented in Section 
[T] (and assumed solar metallicity as well as negligible red- 
dening) . 

The models and assumptions provide us with a predic- 
tion of f\, the spectral energy distribution (SED) received at 
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Figure 10. Difference in the K magnitude (i.e. KI-K2) as a 
function of the K magnitude (average of Kl and K2) for all the 
photometric and astrometric candidates (black dots). The error 
bars on Kl— Kl are given as blue vertical lines, while the hori- 
zontal green line marks zero. The variable candidates from our 
Praesepe cluster candidate list are presented as black circles with 
dots. At the bottom of the figure, the NextGen isochrone masses 
are at the bottom while the DUSTY isochrone masses are above. 
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from the source. We need to 
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Table 8. Seven Praesepe member candidates identified as likely variables, with two doubtful candidates. We give RA, DEC, K, KI — K2 
and comment each object after visual inspection of the UKIDSS GCS DR9 ZY J H K1K2 images. 



RA 


DEC 


K 


Kl - K2 


Comments 


(h m S) 


(° ' ") 


[mag] 


[mag] 




8 


48:45.42 


+22:23:37.2 


14.636+0.012 


-0.211 




8 


49:10.96 


+22:12:13.2 


13.681+0.006 


0.153 




8 


38:07.18 


+21:16:36.8 


14.284+0.006 


0.115 




8 


43:53.41 


+21:01:26.3 


16.918+0.066 


0.229 




8 


36:45.00 


+20:08:45.8 


13.734+0.005 


-0.086 


Faint object observed nearby. Possible companion. 


8 


50:49.83 


+ 19:48:36.5 


12.747+0.002 


-0.270 




8 


41:24.84 


+19:57:27.0 


15.205+0.018 


-0.093 




8 


35:13.24 


+ 17:02:06.5 


13.540+0.004 


-0.190 




8 


41:50.08 


+19:06:18.3 


16.443+0.052 


0.260 


White dots on the object is Kl images (i.e. bad pixels or cosmics). 



convert these spectra to magnitudes in the WFCAM filters 
ZYJHK. Denoting as 5a (A) the (known) total transmis- 
sion function of filter A (including the CCD quantum effi- 
ciency and assuming telescope and instrumental throughput 
is fiat), then the flux measured in the filter is 



/a 



J °° hS A (\)dX 
J™S A {\)d\ ' 



(8) 



The corresponding magnitude ttia in the Johnson pho- 
tometric system is given by 



m A 



-2.51og/A + c A , 



(9) 



where ca is a constant (zero point) that remains to 
be determined in order to put the model-predicted mag- 
nitude onto the Johnson system. We determine this con- 
stant for each of the UKIDSS ZY JHK-bands in the stan- 
dard way, namely by requiring that the spectrum of Vega 
produce a magnitude, ttia, of 0.03 mag in V and zero for 
all colours in all bands. Using the Vega spectrum from 
IColina et al.1 [|l996h . the zero-point constants we obtain for 
the UKIDSS ZY JHK-bands are c z = -22.6587 mag, c Y = 
-23.1075 mag, cj = -23.8080 mag, c H = -24.8365 mag and 
ck = —26.0014 mag. Applying the two equations above to 
the model spectra produces a theoretical isochrone in color- 
magnitude space. 

Unlike most of the open cluster surveys that use one 
magnitude or one colour to estimate masses and effec- 
tive temperatures (T e ff), here we take full advantages of 
the five bands available (i.e. ZYJHK) by determining 
the masses and T e g for each ob ject from its SED as in 
iBoudreault fc Bailer- Jone s (2009) (using a least squares fit 
between the measured SED and the model SED from the 
isochrone) with normalisation at J0. See Fig. QT] for an il- 
lustrative example of the resulting fitting for three cluster 
candidates at three different masses. 

The errors on the derived physical parameters (i.e. on 
mass and T^ff) take into account the errors on the derived 
age of the cluster and the errors on the fitting. 



4 We looked at the minimized sum of the square of the residuals 
as a function of the effective temperature, and there is no general 
trend for a better fit to higher or lower temperature. 



8.2 The cluster luminosity function 

In Fig. [12] we present the LF of Praesepe using the Z and 
J-band photometry of our final list of candidates presented 
in Section[4] Magnitude bins are 0.5 mag wide for both LFs. 
We observe a rise in the LF, a turnover (at Z ~16mag 
and J ~ 15 mag), followed by a decrease to lower magni- 
tudes. The turnover in the J-band LF is at one magni- 
tude brighter than obser ved with the cluster candidates from 
IBoudreault et all l|2010h , based on their NextGen pure pho- 
tometric selection. 

The overall shape of the LF, including the position of 
the turnover, remain unchanged with the binning (binning 
moved by 0.25 mag, bin size increased or decreased by a 
factor of 2). 



8.3 The cluster mass function 

The MF, £(log 10 M), is generally defined as the number 
of stars per cubic parsec in the log arithmic mas s inter- 
val log 10 M to log 10 M + dlog 10 M (|Salpeter|[l955l '). Here, 
we define the MF as the total number of objects in each 

0. 1 log 10 M bin per square degree. Since we do not make any 
corrections for binaries we compute here a system MF. This 
decision was taken considering that the binary fraction we 
have estimated in Section 6 does not cover the whole mass 
range of the MF (and therefore, only a few mass bins could 
be corrected), and also considering that all other MFs used 
in our analysis in Section 8.3 are also system MFs. How- 
ever, a MF with the mass interval 0.45-0.07 Mq corrected 
for binarity is presented at the end of this section. 

We present the MF of Praesepe in Fig. 1131 and the num- 
ber of objects in each mass bin in Tab. [9] We observe that 
the MF has a maximum at ~0.6 Mq and then decreases to 
lower masses. We note also that our determination of the 
MF is unique in the sense that all previous determinations 
of the MF of Praesepe show an increase from 1 to 0.1 Mq 
(Fig|14[l. Previous studies include surveys based on pho tom- 
etry (e.g. IChappelle et al. 1 2005; B oudreault et al1l2010h and 
astro metry (e.g. lHamblv et al.l 1 19951 : iKraus fc Hillenbrand! 
l2007f) . 

One possible explanation is the coverage of our survey, 

1. e. the lack of completeness in the core (i.e. within 1.1° 
from the cluster centre). However, considering that the most 
massive objects should be in the cluster's core, a full cov- 
erage would rather increase the observed slope in our MF 
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Figure 12. Luminosity function of Praesepe using the Z (solid 
line) and J (dotted line) band photometry of our final list of 
ZYJHK member candidates. Magnitude bins are 0.5 mag wide 
for both LFs. 
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Table 9. Number of objects in each mass bin. 
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Figure 11. Fitting of the SED for three cluster candidates in our 
sample. Open circles are the flux measurements from the ZYJHK 
UKIDSS magnitudes, while the filled circles are the least squares 
fits obtained using the isochrones. The values of masses and T a g 
are given in each panel, as well as their name following the IAU 
nomenclature (see Tab.|2]l- 



of Praesepe. To rule out this possibility, we compared the 
MF of Praesepe within 1.25°, and beyond 1.25° (Fig. EH, 
where 1.25° corresp onds to the half- mass radius of Praesepe 
ij Adams et al.|[2002 ). We can see that in both regions of the 
cluster, the MF decreases towards lower masses. 

Incompleteness can also be ruled out as an explanation 
for the discrepancy observed. As it can be seen in Fig. 1131 
the completeness limit and the saturation limit are affecting 
the lowest mass bin and the highest mass bin only. The 



M 
(Mq) 



logio M 

(Mo) 



N 

(pc-3) 



log 



N 



(PC" 3 ) 



0.056 
0.071 
0.089 
0.112 
0.142 
0.178 
0.224 
0.282 
0.356 
0.448 
0.564 
0.710 



-1.25 
-1.15 
-1.05 
-0.95 
-0.85 
-0.75 
-0.65 
-0.55 
-0.45 
-0.35 
-0.25 
-0.15 



2 a 

37±15 
67±7 
75±16 
110±2 
115±4 
144±4 
134±14 
123±7 
131±11 
136±1 
42 a 



-0.163 a 

012+ 0146 
u.uiz_ 221 

270+ ' 045 
U - Z ' U -0.050 

0.485±0.010 
0.504±0.014 
0.602±0.013 

_ +0.043 
u - o,1 -0.048 
534+ 024 

U.OD1_ Q3g 
ri =77+0.004 
u '°' '-0.005 

0.067 a 



a Lower limits. 

incompleteness, alone, cannot explain the overall shape of 
the MF. 

One could question the mass determination or the se- 
lection methods as a possible source of discrepancy. How- 
ever, considerin g that we use the same mass determination 
method as from lBoudreault et all (|2010h . and we use a simi- 
lar filt er set and a similar selection method as in lBaker et al.l 
(20100, such an explanation seems unlikely. 

However, we emphasis that we can reproduce the shape 
of previous determinations of the MF if we remove the se- 
lection based on astrometry. This is presented in Fig. 1161 



5 llBaker et al.ll2010t) combined the UKIDSS DR6 with 2MASS 
and the Sloan Digital Sky Survey. 
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where we show both our determination of the MF of the 
Praesepe, with and without the astrometric selection. We 
can see that the difference between the two MFs goes from 
-0.1 dex at 0.5 M© to +0.3dex at 0.1 Mq. The MF we ob- 
tain for the Praesepe without astr ometric selection is als o 
very similar to the one obtained by Chappclle et al. I (|2005h . 
with an increase of the MF of ~0.2dex from 0.5 down to 
0.1 M e . 

We therefore conclude that our determination of the 
Praesepe MF is different to previous work due to our selec- 
tion based on astrometry and photometry, which includes 
five filters covering the SED from 0.85 to 2.3 /im. Surveys 
based only on photometry would appear to suffer from con- 
tami nation by field M dwarfs and background r ed giants 
(e.g. IChappelle et alj|2005l; iBoudreault et al.ll2010l). 

As for the work of iKraus fc Hillenbran d| (|2007h . the 
overall shape of the MF presents a rise of 0.4 dex from 1 down 
to 0.2 Mq. However, one needs to be careful with the low- 
est bin because this survey is ~70% complete at M5 dwarfs 
(corresponding to 0.15-0.2 M© in Praesepe). 

It could be suggested that a larger spatial coverage 
should in principle detect a larger fraction of low-mass 
stars due to mass segregation. H owever, our s patial cov- 
erage is larger than the survey of I Baker et alJ i|201Ch and 
lHamblv et all (|l995h . and we obtain a larger ratio of the 
number of intermediate-mass stars to lower-mass stars. Mass 
segregation does not appear to be the cause of the discrep- 
ancy in the MF in the low-mass star regime. 

We compare our determination of the MF of Praesepe 
with the Hyades ope n cluster in Fi g. 1171 which has a simi- 
lar age of ~600Myr (|Bouvier et al.ll2008n . We see now that 
both clusters present a decreasing MF to lower masses, un - 
like what was previously reported (|Boudreault et al]|201Cf ). 
However, the Praesepe remains different to the Hyades MF, 
with a decreasing MF of only ~0.2dex from 0.6 down to 
0.1 Mq, compared to ~ldex for the Hyades. 

aresent the MF of a Per (85±5Myr, 
and the Pleiades (125±8Myr, 
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of Praesepe. The a Per and Pleiades MFs are obtained using 
the sames survey and selection procedure as our current 
work on the Praesepe, i.e. using the ZYJHK photometry 
and astrometric data from UKIDSS GCS DR9. We observe 
that the MF of Praesepe in the mass range of 0.072-0.6 Mq 
is more similar to the MF of a Per than the MF of the 
Pleiades. This is surprising as dynamical evolution would 
give a MF with a shape more similar or in between to the 
shape of the MF of the Pleiades and the Hyades, which 
are c lusters with respectively ages of ~125 |Stauffer et al.l 
1998) and ~600Myr. These values are closer to the age of 
Praesepe (r = 590 ^° Mvr; iFossati et al.ll2008ri com pared 
to a Per (85±5 Mvr. iBarrado v Navascues et al.lliooi 'l. 

Finally, we present in Fig. [19] the MF of Praesepe 
from this work corrected for binarity (i.e. between 0.45 and 
0.07 M©). For this, we assume that all binary candidates in 
Fig. [9] have a mass ratio of 1.0 for each object. (With the 
photometry alone, it is not possible to make a reliable esti- 
mation of what would be the binary mass ratio.) The shape 
of the MF obtained is similar to the system MF, i.e. a de- 
crease from 0.6 down to 0.1 Mq. However, we stress that this 
single star MF should be used with caution, considering our 
assumption about the binary mass ratio and the fact that 
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Figure 13. The Praesepe mass function. Error bars are ob- 
tained using the errors on the age of Praesepe. The vertical long- 
dashed line is the lowest completeness limit (at K2 = 18.4 mag, 
M = 0.04 Mq). The vertical dash-dotted line is the highest com- 
pleteness limit (at Z = 20.4 mag, Af = 0.07M©). The vertical 
short-dashed line is the lowest saturation limit (at H = 11.4 mag, 
M = 0.62 Mq). The vertical dotted line is the highest saturation 
limit (at J = 10.9 mag or M = 0.85 Mq). Arrows represent lower 
limits, as we are not complete in these mass bins. 



binary candidates are not known for the whole mass range 
of our system MF. 



9 SUMMARY AND CONCLUSIONS 

In this paper we presented the results of a wide field, near- 
infrared study of the Praesepe cluster using the DR9 of the 
UKIRT Infrared Deep Sky Survey Galactic Clusters Survey. 
We performed an astrometric and photometric selection of 
1,116 cluster candidates out of the 218,141 point sources 
detected towards Praesepe. 

Possible sources of contamination include Galactic disk 
late-type and giant stars and unresolved galaxies. We esti- 
mate a contamination of 11.9% above 0.4 Mq, 9.8% in the 
mass range 0. 15-0.4 M©, and 23.8% below 0.15 M©. 

We investigated the binary frequency of Praesepe using 
the photometry and colours from our cluster candidate s. We 
observ e a binary fraction similar to the simulation of Bate 
(|2012h between 0.07-0.1 M@, ~1.5cr difference in the 0.2- 
0.45 Mq mass interval, and significantly lower by more than 
3cr for the mass rang e 0.1-0.2 Mq. On th e other hand, the 
binary fraction from IPinfield et al.l (|2003h are higher than 
our values and those of Bate! 1 20121 ). However, we note that 
two other works focusing on field l ow-m a ss sta rs have also 
derived binary fractions lower than [Bate] (|2012T l. 

We also studied the variability of the Praesepe candi- 
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Figure 14. Mass function of Praesepe from our present work and 
from pre vious work. We also s how the system Galactic field star 
MF from lChabrier et al.l d2005t) as the dotted curved lines and the 
substcllar limit as a vertical d ashed line. We norm alised all the 
MFs to the log-normal fit of IChabrier et al.l l|2005l ) at ~O.3M 
(logM 0.5). 



dates using the two K-band epochs provided by the GCS. 
We identified seven candidate variables, including one in the 
substellar regime. 

We derived the luminosity function of Praesepe in Z 
and J-band here. We observed that the peak of the J-band 
luminosity f unction is one magnitud e brighter than the one 
reported bv lBoudreault et alj (|2010T ). 

Finally, we determined the mass function of Praesepe, 
which differs from previous studies: while previous MFs 
showed an increase from 0.6 to 0.1 Mq, our MF shows a 
decrease. We looked at the MF of Praesepe at two differ- 
ent regions of the cluster, i.e. within and beyond 1.25°, and 
we observed that both regions show an MF which decreases 
to lower masses. We compared our MF of Praesepe in the 
mass range 0.072-0.6 Mq with the ones of the Hyades, the 
Pleiades and a Per. We concluded that our MF of Praesepe 
is most similar to the MF of a Per although they are respec- 
tively of ~85 and ~600 Myr. Even though of similar age, the 
Praesepe appears different to the Hyades, with a decrease in 
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Figure 15. Mass function of Praesepe for the whole survey area 
compared to the MF of Praesepe within 1.25° of the cluster cen- 
tre, and beyond 1.25° from the cluster centre. The dotted curved 
lines, the vertical dashed line, and the normalization of all the 
MFs at ~0.3M Q are the same as in Fig. ITU 
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Figure 16. Mass function of Praesepe from our present work with 
astrometric and photometric selections (bottom) and without se- 
lection based on proper-motion (PM, top). The dotted curved 
lines, the vertical dashed line, and the normalization of all the 
MFs at ~0.3Mg are the same as in Fie. 1141 We also show in the 
top panel of the figure the difference between the two MF. 



the MF of only ~0.2 dex from 0.6 down to 0.1 Mq, compared 
to ~ldex for the Hyades. 
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Figure 19. Mass function of Praesepe including correction for 
binarity (top) and without correction (this work). The dotted 
curved lines, the vertical dashed line, and the normalization of all 
the MFs at ~O.3M are the same as in Fig. ITU 
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